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Abstract 

Knowledge of the properties of the exchange-correlation functional in 
the form jV xc ([p\], j), where p\(r) = A 3 p(Ar), is important when ex- 
pressing the exchange-correlation energy as a line integral E xc [p) = 
Jq d\ J drjv xc ([p\], j) [3p(r) + r.Vp(r)] (van Leeuwen and Baerends, Phys. 
Rev. A 51, 170 (1995)). With this in mind, it is shown that in the low den- 
sity limit limA-+o / p( r )V 2 jV X c([px], f ) d 3 r < 4tt f p(r) 2 d?r. This inequality is 

violated in the local-density approximation. 
PACS numbers: 31.15.Ew, 71.10.+x,31.25.-v 
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According to density functional theory the electronic groundstate energy Eq can be 
expressed as 

E = T s [p] + U[p] + E xc [p] + Jdr p(r)v(r) (1) 

where v(r) is the external potential, T s [p] is the non-interacting kinetic energy, U[p] is the 
classical electron-electron repulsion energy, and E xc [p\ is the exchange-correlation energy and 
p(r) is the charge- density distribution [|l|-|3| . In a recent paper van Leeuwen and Baerends Q 
discussed energy expressions in density-functional theory using line integrals. The exchange- 
correlation energy, for example, can be expressed as 

em - e xc[Pi } = i ^J*-^ — dr ( 2 ) 

where 7(0) = pi(r) and 7(1) = fair). For practical reasons, it is appealing to choose M 

7(A) = p A (r) = A 3 p(Ar) (3) 
Levy has shown that in the low density limit, A — > 0, 

E c [p x }=XD[p} + ... (4) 
where D[p] is finite and independent of A. Further, since E x [p\] = XE x [p] M it follows that 

\im E xc [px\ = 0. 

A — >0 

With v xc ([p],r) = , the exchange-correlation energy can thus be written as 

E xc [p] = J*d\j d 3 r v xc ([p x ], r) [3A 2 p(Ar) + A 3 r.V Ar p(Ar)] (5) 
= dX J d 3 r jv xc ([p x ], [3p(r) + r.Vp(Ar)] . (6) 

Also see discussion by van Leeuwen and Baerends below Eq. (20) in Ref. 0. 

It is therefore of considerable interest to have information on the behavior of tI'xcQpa]: f ) 
as a function of A. In a recent paper Levy and Perdew [0] showed that the functional 

A[p] = limA- 1 ^ c [p A ] + [/[p] (7) 

A — >U 



is convex, in other words 



f(r) s f_f} P } f(r') d 3 rd 3 r'>0 



5p(r)5p(r') ' 



(8) 



for arbitrary f(r) that satisfies 



J f(r)d 3 r = 0. 



(9) 



In particular, for a finite system or for a system where the charge density vanishes at the 
boundary, /(r) = ^4^, (xj = x,y,z) satisfies Eq. (0). Hence 



dxi ' 

r jgfr) 62A ^ d p^') Mr > > o. 

7 cixj 5p{r)5p{r') dx\ 
For any functional Q[p(r)] invariant under translation, it follows that Q[p(r 
Q[p(r+tR)] where R is an arbitrary vector and t a scalar. As a consequence || 
d 2 Q[p(r+tR)} 



(10) 







dt 



(11) 



t=o 



5p(r)5p(v 



/ 


r<L4[ P ]i 




5p(r) 



R.V[R.Vp(r)] d 3 r. 



Now, since is invariant under translational and R is arbitrary, it follows that 

p(r) d 3 r. (12) 



dp{v) 5 2 A[p] dp(rO ^^3^ = _ /■ [ 5A[p\ 



dxi Sp(r)5p(r') dx\ J [dxf Sp(r) 

The convex nature of A[p] combined with translational invariance therefore imposes, from 
Eqs. ( |10|) and (fl2|) , a constraint on the first functional derivative of A[p] : 

d 2 8A[p] 



dxf 5p(r) 



p(r) d 3 r < 



(13) 



Using the expression 



5E xc [p x [p}} 5p x ([p};r') ^ 



5E xc [px 

5p(r) J Sp x ([p\;r') 5p(r 

= J v xc {\p x [p\;r')X 3 5(r-Xr' 
r N 
A 



v xc ([px; x >■ 



(14) 

(15) 
(16) 



it follows from the definition of A[p] and Eq. (|T3|) that 
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pO) J .2 



1 r 

v H {\p),T) +\imjv xc ([px), 



d 3 r < 0, 



(17) 



where vh 



low density limit. 



sum 

Sp(r) 



J \0gr\ d 3 r. Eq. (|17D imposes a constraint on jv xc ([p\\, v) in the 



Substituting V 2 for Eq. ( |17D is also valid and the constraint can be expressed as: 



1 r 

v h ([p],t) + lim-^ c ([p A ],-) 



p(r)V 2 

-4tt I p(r) 2 d 3 r + limy p(r)V 2 ^. c ([p A ], j) d 3 r < 



(18) 



If we separate the correlation and exchange potentials, and take into account that 
v x([p\], f) = ^([p],r) Eqs. ([L7D and (|l|) may be written as 



P( r ), 2 



1 r 

vh{\p],t) + v x (\p],t) +lim-u c ([p A ],- 



c/ 3 r < 0, 



(19) 



and 



p(r)V 2 



1 r 

v x ([p],r) + hm-t; c ([p A ],- 



<i 3 r < 47r y p(r) 2 d 3 r, 



(20) 



respectively. In the special case of two electrons, ^([p],! - ) = — |f#([p], r). Hence, for this 
special case the constraint on v c alone is 

y p(r)V 2 lim^ c ([p A ],^) d 3 r < 2iv J p(r) 2 d 3 r 

In the local density approximation 



A LDA [p] = -\c\ /p(r)f d 3 r+U[p], 



(21) 



Eq. (p~8|) requires that 



4vr y p(r) 2 d 3 r - ~ |c| y p(r)V 2 p(r)^ d 3 r < 0. 



(22) 



This condition is easily violated. For example, with p(r) = Xm(l + r) m e~ r , with x m 
a normalization constant, the value of |c| for which Eq. (|22|) is violated for can be made 
arbitrarily small by increasing the power m. Eq. (|17]) thus imposes a non-trivial constraint 
on the exchange-correlation potential. 
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